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Parameters for shape, mass and virtual mass are described for the wings and bodies 
of 7 spp. The distribution of wing area and virtual mass follow precise allometric 
relationships, and these parameters can be predicted from the position of the centroid 
of area for each wing. Wing mass distributions show less clear trends due to specific 
effects of the pterostigmas. The different wing shapes between Anisoptera and 
Zygoptera are related to differences in their flight behaviours. No difference in wing 
shape was found between ‘flier’ and ‘percher’ groups of dragonfly. A clear relation¬ 
ship exists between the radii of the centre of mass and the second radius of gyration 
for the odon. bodies. The position of the centre of mass relative to the wing bases is 
related to the flight style and manoeuvrability of these insects. 

INTRODUCTION 

The design of dragonfly wings is the result of many compromises between the 
structural and aerodynamic demands required by each species. The extant Odonata 
consist of three suborders: Anisoptera, Anisozygoptera, and Zygoptera, which as 
their names suggest can be categorised by their relative wing shapes. Anisoptera 
are typified by large, strong and fast fliers able to turn rapidly on the spot; their 
fore- and hindwings are broad and of different shapes. In contrast, Zygoptera are 
characteristically smaller and fly slowly in any direction without turning their body 
(BRODSKY, 1994). They have almost identical fore- and hindwings which narrow 
towards their bases. The extant Anisozygoptera are a minor suborder with only two 
extant species and they have almost identical fore- and hindwings similar to the 
Zygoptera, but a body shape and anal appendages similar to the Anisoptera 
(RUPPELL & HILFERT, 1993); these dragonflies are not considered further here. 


1 Current address: Gatty Marine Laboratory, University of St Andrews, East Sands, St Andrews, 
Fife, KY16 8LB, United Kingdom 
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This study addresses whether differences in the general flight behaviour of 
Anisoptera and Zygoptera can be correlated with wing morphology. 

Various aspects of the functional relevance of wing features have been studied 
for Odonata. HERTEL (1966) described the wing of Aeshna cyanea as having a 
pleated structure near the leading edge, with strong veins top and bottom joined by 
thin membranes. These corrugations are deepest near the wing base and become 
shallower towards the tip and posterior wing margin; the combination of longitudi¬ 
nal veins and the membrane between them form a structure that is strong, light and 
resistant under normal conditions to transverse bending, although it is capable of 
local deformations to absorb sudden impacts at the wing tip (NEWMAN, 1982; 
NEWMAN & WOOTTON, 1986; WOOTTON, 1991). Corrugations give insect 
wings strength with scarcely any weight penalty (REES, 1975a), and the aerody¬ 
namic lift produced by these wings is no worse than that from airfoils with a smooth 
profile (REES, 1975b). Dragonfly wings have a serrated leading edge costa 
(HERTEL, 1966) which promotes the transition from laminar to turbulent flow 
over the wing, resulting in an increase in lift generated by the wing (NEWMAN et 
al., 1977). Small spurs have been noted on the surface of odonatan wings (HERTEL, 
1966; NEWMAN et al., 1977; D’ ANDREA & CARR, 1988, 1989). Newman and 
co-workers showed that the height of these spurs was less than that for admissible 
roughness: ‘‘the amount of roughness which is considered ‘admissible’ in engi¬ 
neering applications is that maximum height of individual roughness elements which 
causes no increase in drag compared with a smooth wall” (SCHLICHTLING, 1968). 
Hence it is unlikely that these spurs have any aerodynamic function, a conclusion 
that has been misquoted by D’ANDREA & CARR (1988,1989). NEWMAN (1982) 
concluded that differences in the construction of the wings probably reflect the 
ways in which they are used in flight rather than their efficiency at performing their 
basic structural functions. 

The wings flip over between the up- and downstrokes, and the supinated wing is 
effectively upside-down during the upstroke. Nonetheless, the wing camber also 
reverses so that it is in the correct sense for lift generation on each half stroke. 
ENNOS (1988) has shown that for Diptera the wing veins diverge posterodistally 
from a twistable leading edge, and aerodynamic loading on such a structure causes 
the wing to twist into camber of the correct sense. Odonatan wing veins conform to 
this model (MAGNAN, 1934; WOOTTON, 1991), and so camber generation may 
also be caused by the aerodynamic loading of the wing. 

The pterostigma on dragonfly wings is a pigmented spot close to the leading 
edge near the wing tip. The pterostigma changes the mass distribution of the wing 
significantly, and effectively moves the centre of mass forward towards the torsion 
axis of the wing (NORBERG, 1972). This reduces inertial wing oscillations and 
flutter, thus raising the critical flight speed above which gliding becomes impossi¬ 
ble. Nonetheless, the centre of mass is still behind the torsion axis for the wing, and 
wing twist at supination is probably caused by inertial forces acting on the wing 
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(NEWMAN, 1982). 

A number of studies have measured first order mass and shape parameters that 
describe dragonflies, such as wing mass, wing area, distribution of wing area, wing 
loading, and moment of wing inertia. Allometries of these parameters have been 
compared to scaling with geometric similarity; departures from geometric similar¬ 
ity with changing size can suggest functional reasons for, or consequences of, trends 
in body and wing shape. Caution should be used when comparing the allometry of 
parameters with respect to body mass, as mean masses for a species can differ 
between seasons (MARDEN, 1989; GRABOW & RUPPELL, 1995). 

MAY (1981) found a relative increase in wing length with increasing body mass 
for Anisoptera, whereas NEWMAN (1982) and GRABOW & RUPPELL (1995) 
report the opposite for Odonata in general. Whatever the precise scaling of wing 
length, all three of these studies show that there is good agreement between wing 
area and geometric scaling; thus larger dragonflies have higher wing loadings (body 
weight / wing area); the wing loading is a measure of how much lift a wing must 
generate per unit area. Presuming that the maximum lift coefficients are similar for 
wings of different dragonflies (WAKELING & ELLINGTON, 1997a), then the 
larger dragonflies must move their wings at higher velocities to generate the re¬ 
quired lift. Indeed, flight speed is expected to be proportional to the square root of 
the wing loading (NORBERG & RAYNER, 1987). Field data for dragonfly flight 
speeds show a general increase with body size and thus wing loading (MAGNAN, 
1934; NEVILLE, 1960; RUPPELL, 1989), and the same is true for butterflies 
(BETTS & WOOTTON, 1988; CHAI & SRYGLEY, 1990; DUDLEY, 1990). It 
should be noted that maximum flight speed is additionally governed by factors 
such as maximum muscle power outputs and so is not directly related to wing 
loading. Damselflies do not show a systematic increase in wing loading with body 
mass because the heavier Calopterygidae have a disproportionately large wing area. 
However, the general wing loadings in damselflies are very much lower than those 
in dragonflies (GRABOW & RUPPELL, 1995). 

There are systematic differences between the shape and mass distributions of 
anisopteran and zygopteran wings. Among the Zygoptera (apart from the 
Calopterygidae), the outer half (by length) of the wing accounts for typically 67% 
of the wing^area; the corresponding figures are 60% for Calopterygidae, and 53% 
and 45% for Anisoptera fore- and hindwings, respectively (GRABOW & RUPPELL, 
1995). Various reasons have been put forward for the shape differences between 
dragonfly and damselfly wings based on flight velocity and wing inertia argu¬ 
ments (WOOTTON, 1991). The flapping speed of any chord wise element of the 
wing increases almost linearly from the wing base to the wing tip. For slow flight 
and hovering the wing base has very low velocity, and so there is a large ratio 
between the velocities of the wing tip and base; at fast speeds, on the other hand, 
the wing bases already have the forward velocity of the dragonfly, and so there is a 
lower ratio between the wing tip and base speeds. The lift produced by each seg- 
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ment of the wings is proportional to the square of the velocity of that segment. 
During slow flight, typical of the Zygoptera, the wing base cannot generate much 
lift because of its low velocity, and so the wing area is concentrated towards the 
wing tip where the velocities are greater. For the fast forward flight of Anisoptera 
the proximal wing region can generate useful lift and so has a larger proportion of 
the wing area; the interial costs of accelerating the wing during each halfstroke are 
also minimised by concentrating the area, and thus mass, towards the wing base. 
However, the inertial power of both dragonflies and damselflies can be totally re¬ 
covered regardless of assumptions about the degree of elastic energy storage in the 
thorax (WAKELING & ELLINGTON, 1997c), and so inertial arguments may not 
be appropriate for extant dragonfly species. Nonetheless, a modern dragonfly fly¬ 
ing with petiolate, zygopteran-like wings would incur larger inertial costs that may 
significantly affect the power requirements for flight. 

Aspect ratio is the ratio of twice the wing length to the mean wing chord, and the 
broad, low aspect ratio, dragonfly wings can be explained by their gliding flight 
(WOOTTON, 1991). Minimising drag during gliding is a compromise between 
minimising profile drag and induced drag from the wings (ENNOS, 1989). At the 
Reynolds numbers at which insects fly, profile drag increases with aspect ratio 
whereas induced drag decreases, and thus there is a particular aspect ratio for which 
drag can be minimised. ENNOS (1989) predicts that for gliding insects of 100 mg 
to 1 g the optimal aspect ratio should be between 5.0 and 7.9, respectively. Ennos 
uses ELLINGTON’S (1984b) formula for predicting profile drag, but this is de¬ 
rived for hovering with large angles of incidence and overestimates the drag on 
gliding dragonfly wings (WAKELING & ELLINGTON, 1997a). Additionally, the 
induced drag from a pair of tandem wings is less than for isolated pairs of wings 
(PRANDTL & TIETJENS, 1957), and so the optimal aspect ratios for minimising 
drag during gliding will be greater than those predicted by ENNOS (1989). 

ELLINGTON (1984a) describes how the use of the moments of area, virtual 
mass and wing mass and their non-dimensional radii can further the understanding 
of the aerodynamic properties of different wing morphologies. For example, the 
second moment of wing mass is equal to the moment of inertia for the wing,while 
in a quasi-steady analysis mean lift and profile power on the wings are propor¬ 
tional to the second and third moments of wing area respectively (WEIS-FOGH, 
1973). Ellington’s parameters have superseded the ‘shape factors’ used by WEIS- 
-FOGH (1973) because the similar sets for describing wing mass and virtual mass 
are also useful for describing wing aerodynamics. There is currently a complete 
dearth of data on these morphological parameters for Odonata, and this study starts 
to redress this lack of information. 

Moments of mass can be measured for the body in a similar way to the wings. 
Many insects beat their wings in a stroke plane that has a fixed inclination from 
their body axis, and a high degree of manoeuvrability is obtained with the centre of 
mass near the wing base and with a low moment of inertia, because the body will 
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respond quickly to desired changes in speed and direction (ELLINGTON, 1984a). 
Indeed, SRYGLEY & DUDLEY (1993) have shown that correlations do exist be¬ 
tween manoeuvrability and the proximity of the centre of mass to the wing bases 
for a range of butterfly species. Zygoptera can vary their stroke plane angles with 
respect to their body axis (PFAU, 1986, WAKELING & ELLINGTON, 1997b), 
and so their turning speed is not limited by the body inertia. 

Dragonflies can be categorised into two groups, ‘fliers’ and ‘perchers’ depend¬ 
ing on their predominant flight behaviour (CORBET, 1983). Perchers observe their 
territories from their perch, whereas fliers do so on the wing. NEWMAN (1982) 
noted that fliers typically have longer abdomens than perchers, which increase the 
moments of body mass and thus damp down any turning moments generated by 
the wings; this may be a compromise between reducing manoeuvrability but in¬ 
creasing body stability during hovering, so the eyes are more effective at detecting 
prey. MAGNAN (1934) reports that the centre of mass lies between the fore- and 
hindwings in the flier Anax parthenope; however, he implicitly states that the non- 
-dimensional hindwing base position is 0.35 body lengths posterior to the front. 
Measurements from scale drawings (ASKEW, 1988) suggest that for A. parthenope 
this non-dimensional hindwing base location is only 0.24, and for the 17 species of 
European Aeshnidae the mean position is 0.23, with the hindwing bases never 
being more than 0.28 body lengths from the front. Magnan’s non-dimensional lo¬ 
cation for the centre of mass, occurring at 0.29 body lengths from the front was 
almost certainly behind both wing bases and would thus agree with Newman’s 
observations for fliers. 

Dragonflies can have among the highest flight muscle mass with respect to body 
mass of any insect (MARDEN, 1989), but the values are less for fliers than for 
perchers (MAY, 1981; MARDEN, 1987). In this study the first and second mo¬ 
ments and radii of body mass are measured, as well as the thoracic muscle mass, 
for a range of dragonfly species, in order to help understand dragonfly body 
morphologies. 


MATERIAL AND METHODS 


The anisopteran species Aeshna cyanea, Libellula depressa, L. quadrimaculata, Orthetrum 
cancellation, Sympetrum sanguineum and S. striolatum , and the zygopteran Calopteryx splendens 
were collected at Granchester Meadows, Milton Country Park, Quy Fen, and Wicken Fen in Cam¬ 
bridgeshire between May and August 1993. They were initially stored in dark specimen boxes, cooled 
with ice to about 4°C; all individuals were analysed for their morphology within 24 hours of capture. 

During analysis the dragonflies were anaesthetised with C0 2 to aid handling them, and the mass m 
of the whole insect was measured on a Sartorius Research R200D balance to the nearest 0.01 mg. 
They were later killed by placing in a freezer for at least one hour. 

Relationships were tested for similarity at the 5% significance level. 

WING MASS PARAMETERS. - The left wings of each dragonfly were cut off at their bases and 
separately weighed for their mass m w to the nearest 0.01 mg. They were then placed on a glass slide on 
a square grid with the long wing axis aligned along one of the sides of the grid. Strips 2-3 mm wide 



40 


J.M. Wakeling 


were cut across the wing, starting from the wing base, resulting in n strips (typically 9-12) being cut 
for each wing. Immediately after each strip was cut it was weighed to the nearest 0.01 mg; each strip 
had mass m r and the distance r. from the wing base to the centre of mass for each strip was taken as 
that to the centre of each strip. 

Non-dimensional wing mass m w equals m>N / m . and the non-dimensional radial distance r equals 
, where R is the wing length. The following parameters were calculated from the distribution of 
mass within the wing according to ELLINGTON (1984a): 


Jkth moment of wing mass 


Non-dimensional radius of kth moment of wing mass 


_ n 

m k = 2 X m . r , 

i=i 



(1) 

( 2 ) 


A small mass loss occurred each time a strip was cut resulting in 2m, being typically 8% less than 
the initial mass. The Ath moment of wing mass was corrected for this mass loss by multiplying it by 
the ratio of initial wing mass to 2m, in accordance with ELLINGTON (1984a). 

WING SHAPE PARAMETERS. - The right wings of each dragonfly were cut off at their bases; 
they were each weighed for their mass m w to the nearest 0.01 mg, and maximum length R was meas¬ 
ured with Vernier calipers to the nearest 0.01 mm. The wings were then photocopied at x 2 enlarge¬ 
ment, and this image was taken into a Macintosh Ilex computer via a Panasonic BL600 video camera 
coupled to a Neotech ImageGrabber24 card. A grid of 51 equally spaced parallel lines was superim¬ 
posed onto the wing image, with the lines perpendicular to the major wing axis. The points of inter¬ 
section of this grid and the wing outline were digitised and represent the comers of 50 adjacent 
parallelograms within the wing shape. Each parallelogram has area s- and mean length C f (the mean 
wing chord for the strip). Due to pixellation errors from the image grabber board, square objects were 
grabbed with one side 0.3% longer than the other; this error was less than that for the other procedures 
in morphology analysis, and so was ignored. 

The following parameters were calculated from the distribution of wing area according to ELLING¬ 
TON (1984a). g is equal to gravitational acceleration, p is air density (1.205 kg m 3 ) and p w is wing 
density (1200 kg m 3 ): 


so 


Wing area 

5 = 2 j 

i*l * * 

(3) 

Aspect ratio 

4 R 2 

AR = ~' 

(4) 

Wing loading 

mg 

Pw = *y* 

(5) 

Non-dimensional mean wing thickness 

A m 

* = ^ 

(6) 

kth moment of wing area 

= 2 ! s, r*. 

(7) 

Non-dimensional radius of Jfcth moment of area 

s 

II 

sjp 

(8) 
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Virtual mass 

__ 50 

v - V c * 

(9) 

Non-dimensional virtual mass 

A _ v(ziR) 2 

2 ptrR}' 

(10) 

kth moment of virtual mass 

v -£*l Wl \ 

(ID 

Non-dimensional radius of kth moment of virtual mass 


(12) 


BODY PARAMETERS. - The mass of the body without the wings, was weighed to the nearest 
0.01 g, and its maximum length L to the tip of the genitalia measured with Vernier calipers to the 
nearest 0.01 mm. The bodies were then frozen for at least one hour to make them rigid; the freezing 
procedures produced no perceivable change in mass from m^, and it was assumed that the position of 
the centre of mass also remained unaffected. Video images of the bodies adjacent to plumb lines were 
taken with the bodies suspended by a pin through the top of the first abdominal segment and then with 
a pin through the forewing bases. These images were later analysed by drawing lines through the pin 
positions, and parallel to the respective plumb lines. The lines intersect at the position of the centre of 
mass, which is at distances / from the tip of the head and / j from the forewing base axis. With the pin 
still through the forewing bases, the period of oscillation of the body swinging as a pendulum was 
measured. An infrared light beam was set up so that it was broken on each half cycle of the oscillation; 
a signal from the beam was displayed on an oscilloscope, and the period t for each oscillation thus 
measured. 

Non-dimensional body length ? equals */£,. The moment of inertia of the body and second radius of 
gyration were calculated according to ELLINGTON (1984a): 


Moment of inertia of body / b » 

A 

Non-dimensional radius of gyration of / b / 2 = 

The head, abdomen and legs were finally removed from the thorax with a scalpel, and the thorax 
was weighed to the nearest 0.01 mg. The thorax was submersed in boiled water for 30 seconds and 
allowed to float in the water for a further five minutes. The softened thorax was then bisected between 
the wing bases, and the muscles removed with a pair of forceps. What remained was mainly thoracic 
cuticle, and this was air dried on paper tissue for 24 hours before a final weighing. The thoracic 
muscle mass m was taken to be the diffemce in mass between the fresh thorax and that with the 

m 

muscle removed. This will slightly over-estimate the muscle mass because it will include the mass of 
other tissue such as tracheae that was also removed, but the error should be small. Non-dimensional 
muscle mass ih m - m m/ m . 


1 i'V ,J 

(13) 

4U 2 


(14) 

\ m b L l l 


RESULTS 


Of the seven species studied, Aeshna cyanea is a flier and the rest are perchers. 
The reason for the imbalance partly reflects the ease at catching stationary perch¬ 
ers from their perch compared to catching fliers on the wing. The damselfly 
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Calopteryx splendens is a percher, as are all damselflies, but for the purposes of 
comparison it is treated separately from the anisopteran perchers Libellula depressa , 
L. quadrimaculata, Orthetrum cancellatum, Sympetrum sanguineum and S . 
striolatum. It should be noted that the Calopterygidae are atypical of damselflies in 
that they have a disproportionately high wing area (GRABOW & RUPPELL, 1995). 

WING PARAMETERS 

The measured wing parameters are given in Table I for the forewings and Table 
II for the hindwings; Table I additionally gives the identification code, ID, and the 
body mass and wing loading for the odonates. The data for the non-dimensional 
radii are plotted in Figures 1-4. It is clear from Figure 1 that there are three classes 
of wing shape for these species. Anisopteran hindwings are broad-based, with their 
centroid of wing area within the proximal half of the wing. Anisopteran forewings 
have their wing area more evenly distributed, with the centroid of area just distal to 
the midpoint of the wing; they also show less variation in shape than the hindwings. 
Zygopteran fore- and hindwings are indistinguishable from each other in both shape 
and mass distribution, and the majority of their area is in the distal half of the wing. 

The relations between the non-dimensional moments of wing area, mass and 
virtual mass have not been explained and have been coined ‘laws of wing shape': 
rules that are obeyed even if the reasons for doing so are unknown (ELLINGTON, 


Table I 

Morphological wing parameters for Odonata forewings 
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1984a). Allometric equations that describe the relationships between r 2 (S) and r^(S) 
with r^S), and r 2 (v) with r,(v) are given by: 

? 2 (S) = 0.919r,(5)^™, (15) 

? 3 (5) = 0.886? | ( 1 S) OM * f (16) 

f 2 (v) = 0.925?, (v) 073 *. (17) 

All these relationships are slightly but significantly different from those given by 
ELLINGTON (1984a) for a diverse range of 20 different insect species from six 
orders. It is no surprise that one particular order, Odonata in this case, has slightly 
different wing shape functions from insects in general. Nonetheless, within the 
Odonata there are universal laws of wing shape, given by equations (15) to (17), 
which can be used to predict higher non-dimensional radii from the first non-di¬ 
mensional radius of area and virtual mass. The allometric equation to describe the 


Table II 

Morphological wing parameters for Odonata hindwings 
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0.573 

1.057 

0416 

0482 

1 90 

00124 

0358 

0449 

SS3 

238 

30 90 

54223 

704 

0456 

0527 

0 578 

1.059 

0425 

0490 

2 11 

00119 

0405 

0499 

SS4 

2.25 

28 80 

457 58 

7 25 

0449 

0524 

0577 

1.047 

0412 

0.485 

208 

00128 

0377 

0475 

SS5 

2 16 

29 73 

500.44 

706 

0452 

0.525 

0577 

1051 

0419 

0487 

182 

00118 

0383 

0471 

SS6 

248 

29 70 

48801 

723 

0447 

0522 

0 S75 

1049 

0410 

0481 

1 97 

00153 

0376 

0471 

SS7 

200 

2897 

465 20 

721 

0442 

0518 

0571 

1053 

0403 

0473 

1 83 

00120 

0383 

0483 

SS8 

235 

29 59 

480.08 

730 

0.450 

0.525 

0578 

1.045 

0.414 

0.487 

1.78 

0.0125 

0368 

0464 

SS9 

2 II 

29 68 











0400 

0484 

CSI 

1 81 

2936 

38195 

903 

0 565 

0613 

0649 

1 159 

0 593 

0628 

1.98 

00135 

0416 

0490 

CS2 

165 

28 77 

387 11 

8 55 

0.553 

0604 

0642 

tits 

0 579 

0618 

176 

00124 

0441 

0 508 

CS3 

1 84 

29 28 

40692 

842 

0.551 

0.603 

0641 

1.119 

0 577 

0615 

1.49 

00112 

0444 

0511 

CS4 

1 70 

29 14 

38053 

8 93 

0.557 

0607 

0645 

1.126 

0 585 

0622 

143 

00125 

0425 

0500 

CSS 

1.90 

29 49 

406 70 

855 

0549 

0601 

0640 

1.115 

0 573 

0612 

1 80 

00133 

0439 

0505 

CS6 

1.79 

29 06 

40689 

830 

0549 

0602 

0641 

1.109 

0 576 

0616 

202 


0423 

0491 

CS7 

1.36 

27.04 

34220 

854 

0550 

0603 

0.642 

1.IU7 

0.579 

0.619 

1 66 

00121 

0420 

0487 



















44 


J.M. Wakeling 


relation between r^v) with r\(S) is 

r,(v) = l.592r | (S) uw . (18) 

A combination of equations (15) to (18) can be used to predict the non-dimen¬ 
sional radii of area and virtual mass for these odonates from the centroid of area, 
Pj(S) with a mean error of 0.5% and never more than 2.9%; the relationships from 
ELLINGTON (1984a) would predict these radii with a mean error of 1.1% and 
never more than 2.7%. Both sets of estimations are within the intraspecific varia¬ 
tion for these radii. Figure 3 shows that the allometric relationship between r,(v) 
and Tj( 5) does not fit the data for C. splendens as well as for the Anisoptera. The 
allometric relationships between r,(v) and r,(S) are significantly different for the 
anisoptera and C. splendens ; the relationship for the Anisoptera is also signifi¬ 
cantly different from equation (18) and that in ELLINGTON (1984a), and is 

r 1 (v)= 1.720r 1 CS) ,7# '. (19) 

Predicting the non-dimensional radii of area and virtual mass for Anisoptera 
using equations (15) to (17) and equation (19) results in a mean error of 0.3% and 
never more than 1.4% error. A more detailed study is required if r j(v) is to be 
predicted accurately from r^S) for Zygoptera. No allometric relationship is given 
between r,(v) and r t (S) for Zygoptera as the data available are only from one 
species, C. splendens. 

The non-dimensional radii for wing mass are much more variable than those for 
area and virtual mass. r 2 (m) shows a positive allometry with r { (m), but there are 
different allometries for each species. The reason for this variability is largely due 
to the effect of the pterostigmas on the mass distribution. The pterostigma has been 
described as an inertial regulator of wing pitch by changing the centre of wing 
mass and the moment of inertia of the wing (NORBERG, 1972). Different 
pterostigmas occur for each species, and their effect on the non-dimensional mo¬ 
ments of mass is thus also different for each species. There is no general allometric 
relationship that will predict r 2 (m) from r^m) for Odonata. 


Figs 1-5. The allometric relationships (indicated by the curves): (1) Radii of the second and third 
moments of wing area, r 2 (S) and r 3 (S), plotted against the position of the centroid ?,(£); - (2) Radius of 
the gyration of the virtual wing mass r 2 (v), plotted against its centroid ?,(v); - (3) Positions of the 
centroids of virtual mass r ,(v), plotted against those for wing area Pj(5); - (4) Radius of the gyration 
of the wing mass r 2 (m), plotted against the position of the centre of mass r,(m) [the curves show the 
allometric relationships for each species]; - (5) Radius of gyration f 2 for the odonate body about the 
wing base, plotted against distance f , from that axis to the centre of mass [the solid line shows the 
relationship obtained by least-squares linear regression]. 



r 2 (S) and f 3 (S) 
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BODY PARAMETERS 


The measured body parameters are given in Table III, and the relationship be¬ 
tween f 2 and f j is shown in Figure 5. 

The centre of mass occurred between the wing bases for all the percher 
anisopterans except L . depressa ; this will maximise manoeuvrability by maximis¬ 
ing the responsiveness of the body to turning torques from the wings. L depressa , 
the broad-bodied chaser, is unusual in that its abdomen is relatively short and wide 
compared to the other perchers, and in the male this abdomen is a very conspicu¬ 
ous powder blue; this broad abdomen is probably used as a sexual signal for the 
male at the expense of a loss of manoeuvrability during flight. The centre of mass 
for the flier A. cyanea occurred behind the wing bases, confirming NEWMAN’s 
(1982) distinction between the two behavioural groups. The damselfly C. splendens 

Table III 

Morphological body parameters 


species 

sex 

ID 

~T- 

mm 

m 

mg 

mb 

mg 

m m 

mg 


l 

h 

h 

4 

mg m 2 

Aeshna cyanea 

M 

AC1 

7432 

677 1 

6653 



0 236 

0 094 

0 224 

18547 


F 

AC2 

72 94 

1085 6 

10358 

369 3 

0340 

0304 

0.137 

0 286 

452.33 


F 

ACT 

71 63 

1162 3 

1115.0 

369.9 

0 318 

0 312 

0 149 

0285 

465.19 


F 

AC4 

7849 

818 2 

775.9 

291 3 

0.356 

0267 

0.112 

0 268 

346 42 


F 

AC5 

75.27 

644.7 

612.6 

2182 

0338 

0.253 

0.126 

0 268 

250.05 


F 

AC6 


567.6 

539 5 

1962 

0.346 





LibeUula depressa 

M 

LD1 

4550 

443 6 

422.7 



0368 

0 161 

0.377 

123.99 


M 

LD2 

4546 

471 0 








Ubellula 

M 

LQ16 

46 99 

373 8 

353 3 



0301 

0079 



quadrimacukm 

M 

LQ18 

47 00 

295.0 

262 1 



0.296 

0055 




M 

LQ29 

47 31 


361 6 



0295 

BmI 

0 257 

53.80 


M 


45.97 

377 9 

3542 

183 5 

0486 

0315 

0 128 

0 261 

5067 


M 

LQ31 

46 76 

408 9 

385 5 

1866 


0.324 

0.129 

0 265 

58 99 

Orthetrum 

M 

OC1 

46 83 

4158 

394 1 



0 294 

0082 

0216 

40 18 

cancellatum 

M 

OCT 

48 43 

4491 

426.9 

243.8 

0 543 

0.306 

0.122 

0.257 

66.18 


M 

OCT 

4828 

443 7 

419.3 

233 3 

0.526 






M 

OC4 

49 31 

443.1 

421.8 

227 1 

0513 

0320 

0 122 

0 258 

65.27 


M 

OC5 

4732 

4150 

387.6 

2109 

0 508 

0291 

0 115 

0243 

50.97 


M 

OC6 

46 39 

382.1 

362.0 

1992 

0 521 

0279 

0.093 

0 232 

42 38 


M 

OC7 

49 48 

3507 

332.8 

1794 

0.512 

0 279 


0 184 

27.60 

Sympetrum 

M 


34 15 

1353 

127.4 

655 

0 484 

0256 


0192 

5.46 

sanguineum 

M 

SSan2 

34 29 

121.9 

1153 

600 

0492 

0278 

0075 

0239 

7.77 



SS«n3 

34 81 

114.0 

109 1 

56 1 

0492 

0.273 

0073 

0 200 

534 



SSao4 

33 88 

125 5 

120 8 

61 3 

0489 

0.279 

0063 

0217 

6 52 



SS«n5 

3364 

1330 

125 9 

643 

0483 

0280 

0 041 

0 173 

434 



SSan6 

3164 

1115 

105.3 

535 

0479 

0.255 

0053 

0205 

445 



SSan7 

34 40 

138 5 

1313 

69 9 

0505 

0 270 

0078 

0253 

9 92 


M 

SSan8 

33 30 

1084 

1028 

499 

0460 

0248 

0 032 

0 197 

439 


M 

SSan9 

36.72 

139.3 

131 6 

68 1 

0489 

0270 

0074 

0259 

1186 


M 

SSanlO 

3426 

1336 

1258 

63.4 

0475 

0282 

0077 

0228 

7 69 

Sympetrum 

M 

SSI 

43 12 

153 9 

1423 



0265 

0 065 

0 241 

1673 

striolatum 

F 

SS2 

42 54 

122 9 

113 6 



0.259 

0 088 

0243 

12.11 


F 

SS3 

43 40 

112 7 

102.7 



0268 


0247 

11.76 


M 

SS4 

41.79 

1080 

99.3 



0.272 

0 113 

0.268 

1252 


F 

SS5 

42.32 

1189 

1084 



0270 

0096 

0.243 

11.44 


M 

SS6 

41 81 

1259 

116.1 



0246 

0041 

0207 

872 


M 

SS7 

43 28 

Him 




0.256 

0.066 


977 


M 

SS8 

43 00 

131.7 

121.5 



0266 

0.035 

0 164 

6.15 


M 

CSI 

4528 

91 0 

819 

297 

0327 

0264 

0 089 

0 279 

13 14 


M 

CS2 

4337 

936 

859 

34 3 

0367 

0 311 

0.126 

0330 

1760 


M 

CS3 

47 16 

1236 

113 4 

42 3 

0342 

0 315 

0.146 

0342 

2942 

■ 

M 

CS4 

4836 

119 1 

1088 

42 0 

0353 

0324 

0.165 

0.339 

29.27 


M 

CSS 


105.6 

987 

39 3 

0372 






M 

CS6 

47 97 

882 

804 

298 

0337 

0.278 

0.133 

0319 

18 86 

■ 

M 

CS7 

4474 

81.8 

75.7 

28.9 

0.353 

0.256 

0.103 

0285 

12 24 
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also had a centre of mass behind its wing bases, but this does not necessarily repre¬ 
sent a loss of manoeuvrability because zygopterans can change direction by alter¬ 
ing the position of their stroke planes without needing to change their body attitude 
(WAKELING & ELLINGTON, 1997b). The least-squares linear regression de¬ 
scribing the relation of f 2 with f ( is given by 

f 2 = 1.154 f, +0.142. (20) 

The non-dimensional muscle mass m m was less for the flier than for the percher 
anisopterans. MAY (1981) notes that this is due to both a lower thoracic mass and 
a smaller thoracic volume fraction of flight muscle in fliers. 

DISCUSSION 

In his study on the morphologies of the wings from six insect orders, ELLINGTON 
(1984a) noted that there were general laws of wing shape and that the radii of the 
moments of area and virtual mass could be predicted from r^S) using allometric 
relationships. The allometric relationships between the radii of wing mass and 
r,(S) were not strong enough to be of predictive value, and so he suggested addi¬ 
tionally determining the position of the centre of mass, r^m), by balancing the 
wing on a knife edge. Data for these radii for six butterfly species (BETTS & 
WOOTTON, 1988) showed allometries that were not significantly different from 
those in the ELLINGTON (1984a) analysis, thus reinforcing the idea of general 
laws of shape. Further studies on butterfly morphologies (DUDLEY, 1990; BUN¬ 
KER, 1993) revealed shape and virtual mass allometries that were significantly 
different from Ellington’s, although the radii of the moments of mass were not. 
Dudley attributed these differences to the fact that his analysis was of butterflies 
with wing areas much more skewed distally than for the insects previously studied. 
Dudley then states that “... while wing mass distributions may be conservative across 
the Insecta, wing area distributions are not necessarily so”. The present study shows 
that the odonatan allometric relationships for the radii of area and virtual mass are 
significantly slightly different from the general relationships given by Ellington, 
but that the radii for the moments of mass are very variable and depend on the 
nature of the pterostigmas in the wings. Whereas the relationships derived from a 
broad range of insects (ELLINGTON, 1984a) are a good general starting point for 
predicting the radii of the moments of morphological parameters, it should be re¬ 
membered that within a particular insect order the wing shapes may differ from the 
norm. Particular care should be exercised when using allomteric relationships to 
predict wing parameters which may be affected by the peculiarities of particular 
insect groups, for example butterfly planforms or dragonfly mass distributions. 

The moment parameters of the wings give a global description of their shape. 
There are times, however, when an analytical description of shape is required, for 
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example during a quasi-steady analysis of flight aerodynamics (WAKELING & 
ELLINGTON, 1997c). ELLINGTON (1984a) describes a Beta distribution that 
can be used to predict the normalised chord c = c f^ in terms of radial position r. 
This distribution can be expressed by 


A 

c = 


r'-'O-r)*-' 


T(p+q) , 
IWfo) 


( 21 ) 


where T() is a Euler gamma function, and the parameters p and q are the following 
functions of the moments of area, 



The Beta distribution is compared with the measured values for c in Figure 6 for 
both anisopteran and zygoptera fore- and hindwings; it can be seen that the distri¬ 
butions give remarkably good fits across the entire range of r regardless of wing 
shape. 


(a) __- 





p.._- 



Fig. 6. Comparison between the non-dimensional chord ( c ) measured along the wing length ( r ) 
(dashed line) and the Beta distribution obtained by matching moments (solid line), for the four types 
of wing: (a) Sympetrum sanguineum forewing (SSanl); - (b) 5. sanguineum hindwing (SSanl); - (c) 
Calopteryx splendens forewing (CS1); - (d) C. splendens hindwing (CS1). - [Abbreviations follow 
those in Tabs MIL] 
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WING MORPHOLOGIES FOR ANISOPTERANS AND ZYGOPTERANS 

The precise phylogeny of the Odonata is still under debate (for example FRASER, 
1954; HENNIG,1981; PFAU, 1991) and it is not the role of this study to enter into 
this discourse. It is generally accepted, however, that the Anisoptera are a 
monophyletic group which is derived from zygopteran lineage. The zygopteran 
wing form of identical fore- and hindwings, also seen in the Anisozygoptera, is the 
ancestral condition, and the Anisoptera have later developed modified wings that 
are no longer identical. Each wing form has its advantages, and these are outlined 
below. 

The zygopteran wings, which stem from a narrow base and enlarge to concen¬ 
trate their area in the outer half, show a good adaptation to the clap and fling 
wingbeats which the damselflies typically perform (RUDOLPH, 1976; MARDEN, 
1987). The clap and fling mechanism is peculiar to zygopterans, and is the reason 
for the high aerodynamic performance that they achieve from their wings 
(WAKELING & ELLINGTON, 1997b, 1997c). To perform a perfect fling the wings 
begin with their surfaces touching and then fling open, leading edge first, until they 
finally come apart as the trailing edges separate. For insects with a small wing base 
separation, like the zygopterans, the wings need to taper into a narrow base if they 
are to maximise the contact time of the trailing edge during the fling; hence the 
petiole nature of their wing planform. 

Anisopterans are generally larger than zygopterans, and their higher wing load¬ 
ing means that they typically fly faster. Frequent flights at higher speeds benefit 
from the wing area being more evenly distributed along the span than is the case 
for zygopterans. The reason for anisopteran forewings having a different planform 
to their hindwings can again be attributed to wingbase separation. The gap be¬ 
tween the fore- and hindwing bases is relatively small, and the forewing is unable 
to have an enlarged wingbase without physically interfering with the hindwing; 
additional wing area cannot be placed ahead of the leading edge because this acts 
as a structural girder which could not support additional frontal area. The hind wing 
is able to increase its area near the base, and this has been done to a greater extent 
than in the forewing. 

WING MORPHOLOGIES FOR FLIERS AND PERCHERS 

To draw conclusions on differences between the wing morphologies for 
anisopteran fliers and perchers requires more information because only one flier 
species, A . cyanea , was caught during the study. CORBET (1983) states that mem¬ 
bers of the Aeshnidae and Corduliidae are typically fliers, and the Gomphidae and 
Libellulidae are typically perchers, and so values of r^S) for seven species from 
each of these families were compared to test whether fliers have different wing 
planforms from perchers. Where r,(S) was not available from species caught in 
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this study, the wing planforms from ASKEW (1988) were used; he traced these 
from projected images to preserve correct shape and proportions. The results are 
shown in Table IV. The different families have characteristic wing planforms (as 
measured by r^S)) but the planforms are not significantly different between fliers 
and perchers. The different flight behaviour of fliers and perchers has not led to the 
evolution of differently shaped wings for each group. 

Table IV 

Non-dimensional radii for the centroids of wing area for Anisoptera; data from this study and draw¬ 
ings in ASKEW (1988) 


Behaviour Family 

Species 

Forewing r ( (S) 

Hindwing r,(5) 

FLIERS Aeshnidae 

Aeshna juncea 

0.51 

0.48 


A. cyanea 

0.51 

0.47 


A. grandis 

0.53 

0.49 


Anax imperator 

0.50 

0.48 


A. par the nope 

0.51 

0.48 


Brachytron pratense 

0.52 

0.49 


Boyeria irene 

0.53 

0.48 

Corduliidae 

Cordulia aenea 

0.52 

0.48 


Somatochlora metallica 

0.51 

0.47 


S. alpestris 

0.50 

0.47 


S. flavomaculata 

0.51 

0.47 


Epitheca bimaculata 

0.49 

0.44 


Oxygastra curtisi 

0.51 

0.46 


Macromia splendens 

0.50 

0.47 

PERCHERS Libellulidae 

Ubellula quadrimaculata 

0.51 

0.45 


L depressa 

0.50 

0.45 


Orthetrum cancellatum 

0.50 

0.44 


Sympetrum striolatum 

0.51 

0.45 


S. sanguineum 

0.51 

0.45 


Leucorrhinia dubia 

0.52 

0.47 


L rubicunda 

0.51 

0.45 

Gomphidae 

Gomphus flavipes 

0.52 

0.48 


G. vulgatissimus 

0.52 

0.49 


G. pulchellus 

0.53 

0.49 


G. graslini 

0.52 

0.49 


Ophiogomphus cecilia 

0.50 

0.48 


Onychogomphus forcipatus 

0.52 

0.48 


O. uncatus 

0.52 

0.49 
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